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ELASTIC SCATTERING OF ZERO ENERGY 
POSITRONS BY HELIUM ATOMS 
S. K. Houston 
and 

B. L. Moiseiwitsch 
ABSTRACT 

Kohn’s variational method is used to calculate the scattering length for 
collisions between zero energy positrons and helium atoms in their ground state. 
The trial function used depends on several exponential parameters and up to 
twenty linear parameters and includes explicitly correlation between the posi- 
tron and the electrons of the atom. The trial function also depends on the ground 
state wave function 0 O of the helium atom and we have attempted calculations 
using six different approximations to % . The calculations all seem to diverge 
except when a Hartree Fork approximation to X P Q is used, although with a three 
term Hylleraas % we get the value -0.592a Q for the scattering length. 

We have also calculated the effective electron number Z e f f which gives the 
rate of annihilation of positrons in helium. An error in a previous paper by 

Houston (1968) is corrected. 
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ELASTIC SCATTERING OF ZERO ENERGY 
POSITRONS BY HELIUM ATOMS 

1, Introduction 

In a recent paper, Houston and Moiseiwitsch (1968) have calculated the 
scattering length for elastic collisions between zero energy positrons and 
helium atoms using Kohn’s variational method. The trial scattering wave 
function used allowed explicitly for correlation between the incident positron 
and the target electrons. It contained up to twenty linear parameters and three 
adjustable exponential parameters. However the representation used for the 
helium atom ground state wave function was the analytic Hartree-Fock wave 
function of Green et al. (1954) and this did not allow explicitly for correlation 
between the atomic electrons. It is the purpose of the present paper to investi- 
gate the use of different helium atom ground state wave functions, which allow 
for interelectronic correlation in varying degrees, in a variational calculation 
of the scattering length. 

Now it is well known that Kohn’s variational method gives an upper bound 
to the scattering length when there are no composite bound states of the inci- 
dent particle and the target atom (Spruch and Rosenberg, 1960), provided that 
the exact target wave function is used. We cannot of course do this for helium, 
but we have shown below that Kohn's variational method in this case should 
still yield a good approximation to the scattering length, which will probably be 
an upper bound, provided that the helium ground state wave function used is 
fairly accurate. 


Calculations of the scattering length for elastic collisions between positrons 
and helium atoms have also been performed by Kestner et al. (1905) using a 
pseudopotential model, by Hashino (1966) using the method of polarized orbitals 
and by Drachman (1968) who has used a first order adiabatic correlation trial 
function in both variational and non-variational calculations. 

In addition to the scattering length calculations , the positron annihilation 
rate in helium has been calculated using the scattering wave functions obtained 
in two of the different cases examined. These values are compared with recent 
theoretical and experimental results. 

2, Variational • lethod 

We want to derive Kohn’s variational principle for the scattering length 
for the case where the helium atom ground state wave function has been derived 
using the variational principle for bound states. 

Let us suppose that </> denotes the exact ground state wave function of the 
target atom, having exact eigenergy E q , and that the trial scattering wave 
function for the entire system has the asymptotic form 

*r - -V 1 W 

1 

for large radial distances r x , of the incident positron from the atomic nucleus. 
a T is an adjustable parameter representing the scattering length. Then Kohn’s 
variational approximation to the scattering length is 
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where H is the total Hamiltonian for the whole system and the integration is 
over all the spatial coordinates of the positron and atomic electrons. The trial 
scattering function X V T will contain a T and probably other adjustable parameters 
c i (i - 1 , .n ) » and the optimum values of these parameters are obtained from 
the n + 1 Kohn variational equations 
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Furthermore, since there are no composite bound states of the system of 
positron and helium atom (Gertler et al, , 1968) we have 


where a is the exact scattering length (Spruch and Rosenberg , 1960). 

Let us now consider the case where the target ground state wave function 
and energy are determined using the Rayleigh-Ritz variational method for 
bound states. 

Suppose V is the exact scattering wave function for zero energy positrons 
off helium atoms. Then the integral 
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where the integration is as in (2) above. NowV has the asymptotic form 
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as r t -* ", and so a small change in y has the asymptotic form 

- ‘I* J_ b(o«/' 0 ) - b'/> 0 

r i 

Also if II 0 is the Hamiltonian for the helium atom, then 


(7) 


E o } *0 Ho * o dr ' (8) 

where the integration is over the spatial coordinates of the atomic electrons 
only. Now the Rayleigh-Ritz theorem tells us that a first order change in % 
produces a second order change in E 0 , that is, since 

8 E„ •• 0 (9) 
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we have 
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where &» E 0 ia of the second order of smallness. 


Now let us consider small variations *n I caused by changes in V . We 


have 
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( 11 ) 


We are supposing that(w + 8W) contains an api>roximation to 0 O which has been 
obtained using the Rayleigh-Ritz theorem, and so we must include the small 
change in the energy. 


Expanding (11) we have 
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The last term Is an addition to the usual Kohn expression which is derived 
assuming that ^ and are exact. 

Now, using Green’s theorem, we have 
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where the integration dS r is over a large sphere in r x space, n t » being a unit 


vector normal to this sphere. 


When we insert the asymptotic forms (6) and (7) of 'P and -V , the right 
hand side of (13) becomes 


/* 
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Now 


(H-E 0 )¥=0 (15) 

since ¥ is an exact solution, and 

f dr' = 1 (16) 
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so we can write 
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If we let 
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( 18 ) 


for * e ¥ ,0 O , e q and a, in turn, we have 


a ~ a_ + 
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where 
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A = I ip Q dr' , 


(20) 


If we neglect the second order terms (to this order A= 1) then we have the 
usual Kohn variational approximation a v as given by (2). 
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The first of the second order terms 



V(H - E 0 ) bVd? 
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is a positive definite (or zero) number and if the other second order term were 
not there, then the neglect of (21) would lead to the upper bound a<a v . But 
the integral in the second term 



(22) 


is infinite and so we may not even get a good approximation to a. However, we 
can show that fo^ large r lt the two second order terms cancel out. 
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where the integration dr, is over the coordinates of the positron and 
dr - d.'jdr . Now we can choose (0 O + £0 O )such that 
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and, using (16), we have 
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0 n 80 n dr r 0 


Hence for large y jthe second order terms in the variational approximation (19) 


cancel out. 


This implies that for most good choices of 0 Q we shall get a good approxi- 
mation to the exact scattering length. 

3. Details of the wave function 

The trial scattering wave functions used were similar to that used by 
Houston and Moiseiwitsch (1968) and were of the form 
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where the suffix i refers to the positron, and 2 and 3 to the atomic electrons. 
ib is given by 
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(28) 


where 

Xi ( r i* r 2 5 r i 2 ) = r f r 2 mr i2 n >o). (29) 

Our calculations were performed for S= 4, 10, and, in several cases, 20, 
corresponding to the inclusion of all terms satisfying i , m, n, < 1, 2, 3, 
respectively. 

0 OT (r 2 ,r 3 ) is the ground state wave function for the helium atom and we 
used the following six approximations: 
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I 1 

II 2 

III 3 


where 


term Hylleraas 
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The details of these six wave functions are given in Table I. 

Also given in Table I are the ground state energies E QT obtained with these 
wave functions, and values of. 


A ~ 



/ 


where 0 OT is the trial ground state wave function and % is the twenty term 
ground state wave function of Hart and Herzberg (1957). This 4> 0 is so much 
more exact than any of the above six that the values of A given in Table I are 
a very good approximation to the A defined by equation (20). 
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4. Scattering Results 

Table II gives the values of the variational approximations to the scattering 
length obtained using the approximations (30) to »/» 0 in conjunction with (27). 
Table III gives the values of the exponential parameters o,a,0 appropriate 
to each case. 

When the very simple Hylleraas ground state wave function I was used, 
the minimum value of 0.061 a Q was obtained for the scattering length with 
S “ 4 for the values of & and a given. The parameter 0 was fixed at 0.0 
because it was found that the calculation was unstable for non-zero values of 
0 , giving large negative values for a y . With the same values of and/3, S 
was increased to 10 and this resulted in the value -8.397a Q for a . This is 
obviously too large in absolute value and so we conclude that wave function I 
is inadequate for this type of variational calculation. It is worth noting that 
the value of A for this wave function is the worst of the six values obtained. 

In the two term Hylleraas wave function, II, we picked the same values of 
b ,a and 0 that were "best” for I and we noted a similar behaviour in the values 
obtained for a v . There is a marked improvement in the value obtained with 
S = 10, but it is still too large and negative even though a considerably better 
value for A is obtained. 

However with the three term Hylleraas wave function, in, we have obtained 
our "best" overall value for the scattering length. The values of $,a and/3 
were those which gave the best value for a y at the S = 4 stage. This was 
-0.186a 0 and when S was increased to 20, keeping the same values for S.aand 
0 , we obtained -0.59280 for the scattering length. However, as can be seen 

12 
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from Table II, it is not obvious that a y is converging for this trial function 
as S is increased and so it may happen that a v will become large and negative 
after a few more terms, as happened for functions I and II. Nevertheless this 
value of -0,5928,5 is in good agreement with those calculated by Drachman 
(1968) who gets -0,659a 0 using a non-variational method and -0.511a 0 using a 
variational bound method, with Kestner et al. (1965) who get -0.575a 0 , and with 
Hashino (1966) who obtained -0.546a o . It is worth noting that the value obtained 
for A with wave function III is the best we have calculated and only differs from 
unity by 10“ \ 

The results for the Hartree Foe k wave function, IV, have been reported 
previously (Houston and Moiseiwitsch, 1968). In this case the values of b ,a 
and/3 are those which are "best" for S= 10. Furthermore as S is increased, 
the values for the scattering length seem to be converging. In the case of the 
Modified Hartree Fock wave function, V, we only calculated one value for a v , 
using the same values for b,a and/3 as were used for wave function IV. Further- 
more we only took S = 10. The value obtained for a v was not as good as the 
corresponding value obtained using the Hartree Fock function and so this choice 
for the exponential parameters was a very poor one. Had we optimized these, we 
should have gotten an improvement over IV. However this calculation took a 
long time and we did not proceed with it. 

Finally the Hylleraas Eckart wave function VI produced values for \ which 
were disappointing in that they seemed to be converging until S was about 17 
when they suddenly became large and negative, in a similar way to I and II, 
although not to such an extent. In this case it was also found that the parameter 
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// could not bo varied from 0.0 and the values for *? and ' are those which were 


’•best" for S - 4. 

So it seems that with the type of scattering wave functions given by (27) , 
the only reliable form for «/. is IV, the Hartree Fock representation. 
Certainly a better value for the scattering length is obtained using the three 
term Ilylleraas wave function in, but it is not clear that the values obtained 
are converging with increasing S as they seem to be doing when wave function 
IV is used. 

5. Positron Annihilation 

The rate of annihilation of positrons in helium gas is given by the express- 
ion 


" r.. 


N Z 


eff 


(Ferell, 1956) where r Q is the classical electron radius e 2 'mo 2 , N a is the 
number density of atoms in the target gas and Z 6{( is an effective election 
number given by 


Z ef f = j dT IVCW^r 


X [Mr t - r 2 ) + 8(r x - r 3 )] 


(31) 


Z eff is a measure of the probability that the positron and an atomic electron 
are at the same point in space. 
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We have used the total scattering wave function (27) with two of the 
different forms of /' 0 to calculate Z r ff . 

With the Hartree Fock function IV we get the value 3,28 for z ({ and with 
the three term I-Iylleraas function III we get 3,21. These values are for the 
scattering wave functions obtained with S « 20, The results are presented in 
Table IV along with the values for Z eff obtained using the wave functions 
with S = 4 and S - 10. Houston (1968) has previously calculated Z ef f using 
the Hartree Fock form of ./* Q but there was an error in his computer program 
and the published value of 4.61 is wrong. 

Drachman (1968) obtains the value 3.66 for Z ef{ for zero energy incident 
positrons and these may be compared with the recent experimental value of 
of 3.628 ± 0.041 obtained by Lee et al. (1969) for positrons with thermal 
energies. Other experimental values are 3.2 (Daniel and Stump, 1959), 3,25 ± 
0.22 (Duff and Heymann, 1962), 2.46 (Osmon, 1965), 3.84 (Tao and Kelly, 1969), 
and 3.680 ± 0,025 (Leung and Paul, 1969). 

We see that our two values are very close to each other and are also in 
good agreement with experiment showing that our scattering wave functions 
are a fairly good approximation to the exact total wave function. 
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TABLE IV 


Effective Electron Number Z eff 




S = 4 


10 


m 


1.86 


2.39 


20 


3.21 


IV 


2.06 


2.48 


3.28 









